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ABSTRACT 

The accuracy of pairing of the anticodon of the 
initiator tRNA (tRNA fMet ) and the initiation codon 
of an mRNA, in the ribosomal P-site, is crucial for 
determining the translational reading frame. 
However, a direct role of any ribosomal element(s) 
in scrutinizing this pairing is unknown. The P-site 
elements, m 2 G966 (methylated by RsmD), m 5 C967 
(methylated by RsmB) and the C-terminal tail of 
the protein S9 lie in the vicinity of tRNA fMet . We 
investigated the role of these elements in initiation 
from various codons, namely, AUG, GUG, UUG, 
CUG, AUA, AUU, AUC and ACG with tRNAggJ 
(tRNA fMet with CAU anticodon); CAC and CAU with 
tRNA™*; UAG with tRNA™*; UAC with tRNA™*; 
and AUC with tRNA^ using in vivo and computa- 
tional methods. Although RsmB deficiency did not 
impact initiation from most codons, RsmD defi- 
ciency increased initiation from AUA, CAC and 
CAU (2- to 3.6-fold). Deletion of the S9 C-terminal 
tail resulted in poorer initiation from UUG, GUG 
and CUG, but in increased initiation from CAC, 
CAU and UAC codons (up to 4-fold). Also, the S9 
tail suppressed initiation with tRNA^ lacking the 
3GC base pairs in the anticodon stem. These obser- 
vations suggest distinctive roles of 966/967 methy- 
lations and the S9 tail in initiation. 



INTRODUCTION 

The accuracy of mRNA decoding is central to the fidelity 
of translation. Structural studies have revealed that during 
elongation, correct minor groove geometry of the base 
pairs formed by the first two nucleotides of the codon 
with the complementary nucleotides of the anticodon in 
the ribosomal A-site is scrutinized by A 1492, A 1493 and 
G530 of the 16S rRNA (1). In contrast, at the step of 
initiation, although it is clear that the correct pairing 
of the initiator tRNA (tRNA^*) anticodon with the 
initiation codon in the ribosomal P-site determines the 
translational reading frame in an mRNA, it is unclear 
how the pairing between the anticodon and the initiation 
codon is scrutinized. For example, tRNA^J (tRNA^ 1 
with CAU anticodon) initiates not only from AUG but 
also from GUG, UUG and CUG, interestingly, by 
allowing wobble/mismatch at the first position of the 
codon (unlike the third position wobble during elong- 
ation). Also, other than the special cases, such as AUU 
in infC (2) and AUA in mitochondria (3), it is unclear 
whether initiation in vivo occurs from codons with the 
third base wobbling/mismatch. Whether any A-site like 
mechanisms play a role in scrutinizing the interaction 
between the anticodon of tRNA^^ and the initiation 
codon of an mRNA in the P-site (at the step of initiation) 
is unknown. 

High-resolution co-crystal structures of the 70S 
ribosome bound with tRNA 0 * 1 ** in the P-site (together 
with tRNA phe bound in the A- and E-, sites) and an 
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Figure 1. P-site model showing interaction of CAU anticodon (C34, A35 and U36) of tRNA™ 61 with AUG initiation codon (Al, U2 and G3). 
Contacts of G966, C967 and the S9 tail are shown. Spheres associated with G966 and C967 indicate methyl groups at their positions 2 and 5, 
respectively. Lysl27 of S9 tail contacts C34 of the anticodon and is in proximity to G966 and C967. The model was generated with PyMOL vl.3rl 
using PDB accession number 2J00 (4). 



mRNA (4,5) yielded insights into the distinctive features 
of the P-site. A methylated nucleoside of the 16S rRNA, 
G966, stacks against the ribose of the first base of the 
anticodon (C34) of the tRNA™ 61 . The residue C967, 
also a methylated nucleoside, lies next to G966 and 
is also in proximity to tRNA™* (Figure 1). The G966 
and C967 methylations are carried out by the specific 
methyltransferases, RsmD and RsmB, respectively (6,7). 
The structural data (4) have also revealed that the tail of 
S9 protein extends into the P-site and contacts tRNA^ 61 
at positions 33 and 34 (Figure 1). The C-terminal tail 
sequence of the S9 protein is phylogenetically conserved 
(8). And, although modification of residue 966 is also 
conserved, neither the identity of the residue nor the 
nature of the modification is conserved (9,10). The 
residues G966 and C967 lie in the 970 loop, which forms 
helix 31, and mutation of either has been reported to result 
in a slight increase (107-127%) in production of a reporter 
protein (11). The C967 methylation is conserved in 
bacteria, and the two methylated nucleosides are also 
reported to contact the S9 tail via the backbone inter- 
actions (12,13), thus forming a network of interactions 
with tRNA 0 ^ 61 in the P-site. In recent years, roles of 
post-transcriptionally modified nucleosides have been 
studied by several groups, and methylated nucleosides 



have been shown to impact fidelity of initiation (14), 
ribosome recycling (15,16), ribosome biogenesis (17) and 
response to nascent peptide (18). 

To investigate the role of methylations of the adjacent 
nucleotides G966 and C967 and the C-terminal tail of the 
ribosomal protein S9 in initiation, we deleted RsmB, 
RsmD and the last three amino acids of S9 (independently 
or together) from the Escherichia coli chromosome. The 
strains so developed were then used to study initiation 
from a variety of initiation codons with the native or 
mutant initiator tRNAs. We show that the S9 tail and 
the modified G966/C967 contribute to the correct 
poising of tRNA^^ in the ribosomal P-site in E. coli. 



MATERIALS AND METHODS 

Strains, plasmids and DNA oligomers 

E. coli strains, plasmids and DNA oligomers used in this 
study are listed in Supplementary Tables S1-S3. Bacteria 
were cultured in Luria-Bertani broth (LB) or LB-agar 
(LB containing 1.8% agar, Difco) at 37°C or as indicated 
with constant shaking at 200 rpm. Ampicillin (Amp, 
100 ugml -1 ), kanamycin (Kan, 25 ugml -1 ) or tetracycline 
(Tet, 7.5 ugml -1 ) were used as needed. 
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Generation of knockout and knockin strains 

Details of the strains generated are provided in the 
Supplementary Material. 

Generation of chloramphenicol acetyltransferase reporter 
constructs containing mutant initiation codons and the 
corresponding anticodon sequences in the initiator tRNA 
(encoded by metY) 

Reporter constructs were derived from pCAT AUG m^/F 
or pCATuAG^^cuA (19) using two steps of site- 
directed mutagenesis to first mutate the initiation codon 
of the chloramphenicol acetyltransferase (CAT) reporter 
and then the anticodon sequence in the metY gene. 
Reactions (50 ul) contained MOOng DNA template, 
200 uM dNTPs, 20pmol of the required forward and 
reverse primers (Supplementary Table S3) and 1 U 
Phusion M DNA polymerase (2Uul _1 , Finnzymes). All 
constructs were verified by DNA sequencing. 

Preparation of the cell-free extracts 

E. coli strains were grown (3-5 replicates) in 3 ml of LB 
with Amp at 37°C to an OD 600 of -0.6-0.7. The cells were 
harvested, and the cell-free extracts were prepared, 
quantified and stored at -20°C (14). 

CAT assays 

CAT assays using the cell-free extracts (0.3-2.0 jig total 
protein) were performed in 30-ul reactions containing 
500 mM Tris-HCl, pH 8.0, 150uM chloramphenicol, 
0.025 |iCi (430pmol) [ 14 C]-Cm (specific activity 57.8 m 
Cimmol -1 ; Perkin Elmer) and 432 uM acetyl-CoA at 
37°C for 20min followed by separation of the substrate 
and products by thin layer chromatography (15). The 
pixel counts were estimated for the spots corresponding 
to the substrate (S, chloramphenicol) and products 
(P, 1 -acetyl and 3-acetyl chloramphenicol) using the 
Multi Gauge software. Per cent CAT activities were 
calculated as [P/(S + P)] x 100. This value was divided by 
the amount of total protein (in micrograms) used for each 
assay. The activity so obtained was then normalized to the 
relative value of [3-lactamase activity of the corresponding 
extract to avoid any variations in the CAT values because 
of plasmid copy number changes. 

P-Lactamase assay 

(3-Lactamase assays in 125-jil volumes contained 
CENT A™ (0.2 mM, Calbiochem), a chromogenic sub- 
strate, and cell-free extracts (0.8-1.0 ug total proteins) 
for 30min in 10 mM Tris-HCl (pH 8.0). The reaction 
was stopped by addition of 3 ul of 5% sodium dodecyl 
sulphate, and the volume was made up to 400 ul with 
water. Absorbance was measured at 405 nm. The value 
of E. coli SA extracts in each replicate was considered as 
unit, and the values of other extracts were calculated 
relative to this value. 

Calculation of fold differences in CAT activities 

CAT activities of 3-5 replicates were determined. Fold dif- 
ferences for the CAT values of all strains (i.e. E. coli SA and 



its derivatives) for each of the CAT constructs were 
calculated by dividing them with the corresponding value 
of the SA strain for each of the replicates. The averages of 
the fold differences for each strain were calculated and 
plotted along with SEM. Values for SA in such calculations 
remain as unit for each construct. The graphs plotted 
represent CAT activity per micrograms (normalized to its 
relative P-lactamase activity) of each strain relative to the 
control E. coli SA parent strain. A schematic of calculation 
of fold difference in CAT activity is given in Supplementary 
Figure SI. In all assays, a difference of >2-fold (with 
respect to the control) was considered significant. 

Molecular dynamics simulations 

Four codomanticodon pairs, namely, the wild-type AUG: 
CAU, and an array of the mutants CACGUG, CAU: 
GUG and AUC:GAU were used for modelling and 
molecular dynamics (MD) simulations in the presence or 
absence of the S9 tail (the three C-terminal Serl26, Lysl27 
and Argl28, SKR, residues; the S9 tail deletion set is the- 
oretically equivalent to the S9 protein being absent from the 
MD simulations). The region representing the mRNA: 
tRNA 0 ^ complex in the presence of the S9 tail was 
carved out of the crystal structure of the ribosome (PDB 
accession number 2J00) from Thermus thermophilus (4). 
The corresponding structure of E. coli ribosome was not 
chosen because the E. coli structure is complexed with 
tRNA phe instead of tRNA™* (5). However, the region of 
our interest (i.e. P-site tRNA + mRNA + S9) shows high- 
structural similarity to T. thermophilus [Root Mean Square 
Deviation (RMSD): 1.05 A]. The mutant pairs were 
generated using the 'mutate_bases' module in the 3 DNA 
package (20), followed by minimization, equilibration and 
MD simulations for 10 ns each. MD simulations were 
carried out in aqueous medium using TIP3P water model 
using the force field ff99.rna. Na + ions were added to each 
system to neutralize the net charges using tleap in 
AMBER9 (21). The solvation box is 14 A from the 
farthest atom along any axis. The simulations were per- 
formed under normal pressure and temperature (NPT) con- 
ditions with the Berendsen temperature coupling scheme 
and periodic boundary conditions. Particle Mesh Ewald 
(PME) summation was used for long-range electrostatics. 
The van der Waals' cut-off used is 10 A, and the pressure 
and temperature relaxations were set to 0.5ps _1 . A time 
step of 2 fs was used with the integration algorithm, and 
the structures were stored every 1 ps. The trajectories from 
1-10 ns were used for the hydrogen bond analysis. The 
hydrogen bonding analysis was done on each snapshot 
[distance cut-off for identifying hydrogen bonds 
(H-bonds) being 3.5 A]. The H-bonds, which appear in 
>40% of the MD ensemble, were termed dynamically 
stable and were chosen for further analyses and inference. 

RESULTS 

In vivo system for initiation from different 
codon:anticodon pairs 

The in vivo assay is based on a plasmid containing a CAT 
reporter gene (19) whose initiation codon (AUG) was 
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mutated to UAG, CAC, CAU, UAC, GUG, UUG, CUG, 
AUU, AUA, AUC or ACG. The initiator tRNA 
(tRNA™ 61 ) gene (metY) with a corresponding change 
in its anticodon (CAU) to CUA to decode UAG; GUG 
to decode CAC and CAU; GUA to decode UAC and 
GAU to decode AUC was also present on the same 
plasmid. The AUG, GUG, UUG, CUG, AUU, AUA, 
AUC and ACG codons were decoded by the native 
tRNA^u (Supplementary Tables S4 and S5). The initi- 
ation activities of the CAT reporters were monitored by 
growth of the strains on chloramphenicol (Cm) plates, 
CAT assays or immunoblotting (22,23). 

As anticodon (CAU) is the major determinant for 
aminoacylation of tRNA™* by MetRS (24,25), muta- 
tions in anticodon (needed for initiation with many of 
the non-AUG codons) often result in aminoacylation of 
the mutant tRNA™ 61 with non-methionine amino acids 
(Supplementary Table S6). Although, in E. coli, initiation 
with a variety of non-methionine amino acids has been 
reported, their initiation efficiencies vary (26). Thus, to 
avoid any variations because of initiating amino acids, 
the CAT activities were used only to determine the 
ratios (fold differences) of initiation between the strain 
deficient for the specific ribosomal P-site feature(s) 



versus the control strain (wild-type for the same ribosomal 
P-site features) independently for each of the codomanti- 
codon pair. It may also be said that as in our assay design, 
only the initiation codon of the reporter and/or the anti- 
codon of the initiator tRNA are changed, the assay 
reports specifically on the impact of the ribosomal P-site 
elements on initiation. 

Use of initiator tRNAs with GUG or QUG anticodons 

The native tRNA™* with CAU anticodon possesses 
cytosine, a pyrimidine, at position 34. To address 
whether a purine at this position, such as in GUG or 
QUG codons [where Q represents queuosine, a bulky 
modification of guanosine found in tRNAs with GUN 
anticodon, (27)], could be allowed, we introduced the 
reporter pCAT C AC^<^GUG into E. coli BW and its 
Atgt::kan derivative referred to as E. coli SA or SA 
(deficient in Tgt, which is required for Q modification). 
E. coli BW and E. coli SA served to study initiation with 
tRNAQuc and tRNA GUG> respectively. The tRNAs 
encoded by metY G u G are novel variants that initiate, 
but the amino acid attached to these has not been 
identified. As shown in Supplementary Figures S2-S4, 
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Figure 2. Initiation from CAC and CAU codons by tRNAj^Q (encoded by metY GVG ) as assessed by growth on Cm plate (A). Saturated cultures of 
E. coli SA and its ArsmB, ArsmD, and ArsmB ArsmD derivatives harbouring pCAT CAU me/ F GUG (sectors 1-4) or pCAT CAC me/ F GUG (sectors 5-8) 
were streaked on Amp 100 (left) or Amp 100 Cm 250 (right) and incubated at 37°C for ~15h. CAT assays (B) to assess initiation activities of CAC:GUG 
(panel i), CAU: GUG (panel ii), AUG: CAU (panel iii) and UAG: CUA (panel iv) codon: anticodon pairs. Star symbols in the boxes indicate G966 
(anticodon proximal) and C967. Fold differences in the CAT activities with respect to SA (C, taken as 1) in its ArsmB, ArsmD, and ArsmB ArsmD 
derivatives are shown. For reference, the average CAT activities in SA (C) for the CAC: GUG, CAU: GUG, AUG: CAU and UAG: CUA codomanti- 
codon pairs were 1462 ± 374, 493 ± 64, 7314 ± 923 and 5246 ± 763pmol Cm acetylated per 1 jag of total cell-free extract in 20min at 37°C, 
respectively. 
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although the metY G u G encoded tRNA et was queuo- 
sinylated in E. coli BW and it participated in initiation, 
to avoid complications because of its partial modification, 
E. coli SA was used to investigate initiation with 
tRNAgftk and all other tRNAs. 

Methylated G966 and C967 suppress initiation 
with tRNAg^ 

The pCATcac^^Jgug and pCAT CA u^^ 7 G ug con- 
structs were introduced into E. coli SA and its derivatives 
to measure initiation from CAC:GUG and CAU:GUG 
codomanticodon pairs, respectively. Transformants were 
grown in LB with Amp and streaked on LB-agar contain- 
ing 100 |igml _1 Amp (Amp 100 ) as control and on LB-agar 
Amp 100 with varying concentrations of Cm (Amp 100 
Cm too ~ 250 ). Representative plates with Amp 100 or 
Amp 100 Cm 250 (Figure 2A) show that the pCAT CA u 
metY G uQ supported the growth of the strains deleted for 
methyltransferases in singles (ArsmB or ArsmD) or double 
(ArsmB ArsmD) (sectors 2-4) but not of the parent (sector 
1), suggesting better initiation by tRNA^^ in the absence 
of methylations at positions 966 and/or 967. In the plate 
assay, the transformants harbouring pCAT CAC m^ F GUG 
(sectors 5-8) did not show a perceivable difference 
in growth (because of higher level of initiation from 



the matched CAC codon). However, as revealed by the 
CAT assays (Figure 2B, panels i and ii), deletion of 
methyltransferases in singles or double resulted in 
increased initiation by tRNAj?^ (~2- to 3.6-fold) over 
the SA parent (C) from both CAC and CAU codons 
(see also Supplementary Figure S6). Interestingly, as 
shown in Figure 2B (panels iii and iv), the deficiency of 
these methyltransferases showed no significant changes 
in initiation with AUG (decoded by native tRNA^J) 
or UAG (decoded by near native tRNA^). 

Methylated G966 suppresses initiation from AUA codon 

Deletion of rsmD resulted in increased initiation with the 
CAU:GUG pair (U:G mismatch at the third position 
of codon) to a level (> 3-fold) that was higher than that 
seen for the perfectly matched pair CAC: GUG (~2-fold) 
(Figure 2B). This led us to investigate whether G966 
methylation discriminated against the third position (of 
the initiation codon) wobble/mismatch. To address this, 
initiation from various canonical codons AUG, GUG, 
UUG (28) and CUG (29) and non-canonical codons 
AUU (2), AUA (30), AUC and ACG with the native 
initiator tRNA (Supplementary Table S4) was 
investigated. As seen in Figure 3, on deletion of rsmD, 
initiation from AUA (third base mismatch with 



2.5n 




AUG GUG UUG CUG AUU AUA AUC ACG 



Figure 3. Initiation activities of the native tRNA^j from various codons. Star symbols in the box indicate G966 (anticodon proximal) and C967. 
Fold differences in the CAT activities with respect to SA (C, taken as 1) in its ArsmB, ArsmD, and ArsmB ArsmD derivatives are shown. For 
reference, the average CAT activities in the SA control (C) for the AUG, GUG, UUG, CUG, AUA, AUU, AUC and ACG constructs were 
7314 ± 923, 3629 ± 492, 4019 ± 1000, 906 ± 132, 238 ± 33, 301 ± 32, 352.5 ± 38 and 48 ± 15pmol Cm acetylated per 1 |ig of total cell-free extract 
in 20min at 37°C, respectively. 



4968 Nucleic Acids Research, 2013, Vol. 41, No. 9 



tRNA^j) was at an advantage (~2-fold), but initiation 
from the remaining codons was not impacted significantly. 
Also, deletion of rsmB (alone or in combination with 
rsmD) did not have much impact. 

Generation of the S9 C-terminal tail deletion strains 

The C-terminal sequence (Serl26, K127 and R128, or SKR) 
of the ribosomal protein S9 is conserved, and the last two of 
these contact positions 33 and 34 of tRNA^ 1 (4). Lysl27 is 
also in proximity to the methylated nucleosides G966 and 
C967 (Figure 1). This prompted us to investigate whether 
the C-terminal tail of the S9 protein contributes to 
decoding/fidelity at the step of initiation. A strain, S9A3, 
was generated wherein the C-terminal SKR sequence of S9 
was deleted to disrupt its contacts with tRNA 0 ^ 
(Supplementary Material). As shown in Supplementary 



Figure S7A, although the S9A3 strain and its derivatives, 
ArsmD S9A3, ArsmB S9A3, and ArsmB ArsmD S9 A3, are 
cold sensitive, they grow normally at 37°C. 

S9 tail suppresses initiation from unnatural 
codon:anticodon pairs 

E. coli SA and its derivatives were transformed with 
pCATcac^Fgug, pCATcau^Fgug, pCAT AU G or 
pCAT UAG m^ Yqxja and used to monitor initiation 
(Figure 4). Initiation from the unnatural codons, CAC 
and CAU, using tRNAj?^ showed ~3- to 4-fold increase 
in initiation in S9 A3 and S9A3 ArsmD strains (panels i and 
ii). However, initiation from AUG using tRNAjJ^ or from 
UAG using a near natural tRNA^ e A (panels v and vi) 
showed no significant changes. Combining rsmD deletion 
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(iv) AUC:GAU (codon:anticodon) (v) AUG:CAU (codon:anticodon) (vi) UAG:CUA(codon:anticodon) 




Figure 4. Initiation with CAQGUG (panel i), CAU:GUG (panel ii), UAC:GUA (panel iii), AUC:GAU (panel iv), AUG:CAU (panel v) and 
UAG:CUA (panel vi) codon:anticodon pairs. A tail ending in SKR indicates ribosomal protein S9. Fold differences in the CAT activities with 
respect to SA (C, taken as 1) in its ArsmD, S9A3, and S9A3 ArsmD derivatives are shown. For reference, the average CAT activities in SA (C) for 
the CAQGUG, CAU:GUG, UAC:GUA, AUGGAU, AUG:CAU and UAG:CUA constructs were 1236 ± 151, 590 ± 103, 241 ± 26, 12266 ± 2866, 
7838 ± 610 and 5454 ± 682pmol Cm acetylated per 1 |ig of total cell-free extract in 20min at 37°C, respectively. 



Nucleic Acids Research, 2013, Vol. 41, No. 9 4969 



with S9A3 (ArsmD S9A3) did not increase initiation any 
significantly (Figure 4). 

S9 tail deletion increases initiation with UAC:GUA but 
not with AUCrGAU codonranticodon pairs 

To investigate whether the S9 tail disfavoured other 
unnatural codon:anticodon pairs, initiation from 
UAC by tRNAgf&and AUC by tRNAgJfo (26) was 
examined. Although in our assays, a steady-state accumu- 
lation of the formylated form of tRNA^^ and both the 
formylated and aminoacylated forms of tRNAj?^ 
were undetectable (Supplementary Figure S8), both 
of these tRNAs initiated, suggesting that the respective 
forms must be used as soon as they are generated. 
Interestingly, the CAT assays revealed that deletion of 
the S9 tail resulted in a subtle increase in initiation from 
UAC (~ 1.5-fold); and a combined deletion of the S9 
tail and rsmD resulted in ~2-fold increase in initiation 
(Figure 4, panel iii), suggesting that S9 tail discriminates 
against the codomanticodon pairs unnatural to initiation. 
An initiation assay was also done with AUC using 
tRNAjJ^j, a perfectly matched pair (panel iv). Here, 
deletion of S9 tail did not change initiation (~0. 8-fold, 
Figure 4). This seemed to be an exception (see later 



in the text for analysis of dynamically stable H-bonds). 
In addition, this codomanticodon pair (AUQGAU) also 
suffered a problem because of background initiation from 
AUC using the native tRNAg^ (Figure 3). 

S9 tail deletion decreases initiation from the 
less frequently used alternate initiation codons 
(UUG and GUG) 

Alternate initiation codons GUG, UUG and CUG with a 
mismatch at the first codon position are not uncommon in 
mRNAs and are decoded by the native tRNA^J. CAT 
reporters that initiate from these codons showed that 
although the S9 tail deletion did not impact initiation 
from AUG, it led to a decrease in initiation from GUG, 
UUG and CUG (Figure 5). In the absence of the S9 tail, 
the mismatched codomanticodon pairs may be less stable, 
and the tail may be required to maintain interactions 
between the anticodon and the codon with first base 
wobble/mismatch. Initiation from ACG; and AUA, 
AUC, and AUU with the native tRNA^ (mismatches 
at the second or the third positions of the codon) showed 
that these activities remained unchanged in the S9A3 
strain (Figure 5). 
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Figure 5. Initiation in absence of the S9 tail alone or in combination with deletion of rsmD with various initiation codons using the native initiator 
tRNA (tRNA^u). A tail ending with SKR in the box on the top right indicates ribosomal protein S9. Fold differences in the CAT activities with 
respect to SA (C, taken as 1) in its ArsmD, S9A3, and S9A3 ArsmD derivatives are shown. For reference, the average CAT activities in SA strain (C, 
taken as 1) for the AUG, GUG, UUG, CUG, AUA, AUU, AUC and ACG constructs were 7838 ± 610, 3198 ± 378, 3788 ± 334, 327 ± 15, 
313 ± 19, 346 ± 34, 631 ± 71.6 and 78 ± 9pmol Cm acetylated per 1 [ig of total cell-free extract in 20min at 37°C, respectively. 
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Figure 6. Effect of deletion of the S9 tail on initiation with the anticodon stem mutants of tRNA^f^. Various mutations in the anticodon stem are 
as shown (A). Fold differences in the CAT activities with respect to SA (C, taken as 1) in its S9A3 derivatives are shown (B) for 1GC (AU, panel i), 
1GC (GU, panel ii), 2GC (AU/GU, panel iii) and the 3GC mutant (panel iv). For reference, the average CAT activities in SA for the UAG:CUA/ 
AU, UAG:CUA/GU, UAG:CUA/AUGU and UAG:CUA/3GC constructs were 2970 ± 900, 888 ± 145, 32 ± 6, 238 ± 33 and 88 ± 19pmol Cm 
acetylated per 1 |ig of total cell-free extract in 20min at 37°C, respectively. 



S9 tail deletion affords better initiation with the 2GC and 
3GC mutants of tRNA™ 6 ' 

The S9 tail is flexible, and we hypothesized that it may sense 
alterations in the conformation of the anticodon arm of the 
initiator tRNA. To investigate whether the S9 tail imposes 
any preference for the tRNA^ 61 via its highly conserved 
feature of the three consecutive GC (3GC) base pairs in the 
anticodon stem, we used tRNA 0 * 1 ** mutants (Figure 6A) 
possessing changes at the first, first and third, third or 
all 3GC base pairs (15,31). The CAT reporters were 
introduced into E. coli SA and its S9A3 derivative. As 
expected from our earlier observation (Figure 4, panel vi), 
the initiation activity of tRNAj?^ (with 3GC base pairs 
intact) from UAG did not change on deletion of the S9 tail. 
However, the initiation activities of the 2GC and the 3GC 
variants (A:U/G:U and 3GC; Figure 6B, panels iii and iv) 
showed an increase of > 2-fold. The increases in initiation by 
the 1GC mutants, not unexpectedly, were subtle (Figure 6, 
panels i and ii). 

MD simulations 

The MD simulations of AUG:CAU, CAC:GUG and 
CAU:GUG codomanticodon pairs with or without the 



S9 tail did not exhibit any drastic conformational 
changes at the backbone level of tRNA^ 61 . However, a 
rigorous H-bond analysis along the trajectory provided 
excellent means of capturing subtle conformational 
re-orientations, taking into account the minor rewiring 
at the side-chain level. The AUG:CAU pair exhibits a 
dynamically stable H-bond (indicated by blue line, 
Figure 7, panel i) with the S9 tail, which exists in >70% 
of the snapshots. Another H-bond (indicated by yellow 
line) was observed in fewer (~57%) snapshots. Of these, 
the one indicated in blue was also reported in crystal struc- 
ture of the initiator tRNA ribosome complex (4). The 
deletion of the S9 tail (leading to loss of these weak 
H-bonds) is expected to show no major effects on the 
interaction of native anticodon (CAU) with the AUG ini- 
tiation codon. In the CAC:GUG codomanticodon pair, 
substitution of C34 anticodon position by G34 affects 
subtle conformational variations during dynamics, which 
induced enhanced H-bonding of the tRNA with the S9 tail 
(Figure 7, panel ii). These H-bonds arrest the intrinsic 
dynamics and have a freezing effect on the tRNA, which 
is released on the S9 tail deletion. Thus, the S9 tail deletion 
is proposed to allow the tRNA regain the conformational 
dynamics important for efficient recognition of the CAC 
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initiation codon. As shown in Figure 7, panel iii, for the 
CAU:GUG codon:anticodon pair, along with the 
mutation at the 34 position of the anticodon, a GU 
wobble (at the third position of the codon) is also 



incorporated. Here, we observe not only the formation 
of additional H-bonds (i.e. enhanced tRNA:S9 inter- 
action) but also loss of a non- Watson and Crick H-bond 
between the third G:U pair in the presence of the S9 tail. 



(i) AUGrCAU 




Figure 7. Pictorial representation of the dynamically stable H-bonds on the tRNA-mRNA and S9 complex from PDB accession number 2J00 in 
AUG:CAU (panel i), CAC:GUG (panel ii), CAU:GUG (panel iii) and AUC:GAU (panel iv) codomanticodon pairs in the presence and absence of 
the S9 tail, respectively. The H-bonds present in 70% or more snapshots are depicted by blue line, whereas those in 60-70%, 50-60% and 40-50% of 
the snapshots are depicted by cyan, yellow and red lines, respectively. The tRNA 0 ^ 61 and mRNA backbones are depicted as deep green wires, and the 
residues of interest (both amino acids and nucleotides), participating in the H-bonds, are shown in stick representation. 

(continued) 
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Figure 7. (Continued). 



On deletion of the S9 tail, the lost H-bond between the 
G:U as also the dynamical properties of the system are 
regained, thus facilitating codomanticodon pairing. These 
observations are in good agreement with the in vivo data 
showing an increase in initiation from CAC and CAU on 
S9 tail deletion (Figure 4). 

Subsequently, as a control, we carried out MD simula- 
tions on the AUC:GAU codon:anticodon pair, which 



showed negligible change in initiation on S9 tail deletion 
(Figure 4). Interestingly, in this pair, S9 tail contacts 
position 33 and 34 of the tRNA (Figure 7, panel iv). The 
same contacts are seen for the wild-type tRNA^^ with the 
S9 tail. Thus, on deletion of the S9 tail, these contacts are 
lost without affecting the conformation of anticodon loop 
of tRNA or codomanticodon interaction. It is, therefore, 
not surprising that initiation with the AUC:GAU pair is 
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not affected on deletion of the S9 tail. Also, the gain of the 
hydrogen bond at position 1 of the codon may compensate 
for the loss of the hydrogen bond made by the S9 tail with 
the positions 33/34 of the tRNA. 

DISCUSSION 

This study reveals the in vivo relevance of the methylated 
G966 and C967 and the C-terminal tail of the ribosomal 
protein S9 in the fidelity of decoding at the step of initi- 
ation. In a specific hypothesis, we propose that, to ensure 
fidelity of initiation, these elements (one or more) allow 
initiation from the authentic initiation codon (AUG), 
facilitate initiation from naturally used alternate initiation 
codons (such as GUG, UUG, CUG and so forth) and 
suppress initiation from unnatural initiation codons 
(such as CAC, CAU and so forth). Thus, in our 
approach, deletion of one or more of these elements 
(methylations at 966 or 967 positions and the S9 tail) 
would result in a minimal change in initiation from 
AUG, a decrease in initiation from the naturally used 
alternate initiation codons and an increase in initiation 
from the unnatural codons. 

As seen for nearly all codomanticodon pairs tested, 
initiation in the strain with a double deletion of the S9 
tail and rsmD was not too different from that deleted 
for the S9 tail alone. Also, in most cases, fold changes 
in initiation were higher for a strain lacking S9 tail than 
the one lacking rsmD, suggesting a more pronounced 
effect of the S9 tail over the methylation at 966. We sum- 
marize our findings in Tables 1 and 2, and we discuss them 
as follows. 

Methylated nucleosides G966 and C967 fine-tune initiation 

A recent study showed that G966 and C967 methylations 
are important at an early step of initiation complex for- 
mation, and their lack led to ~2-fold decrease in binding 
of tRNA™* (32). The importance of G966 in stabilizing 



elongator tRNAs at P-site has also been highlighted (33). 
In our study of their contributions in decoding at the 
P-site during initiation, we observed that the codomanti- 
codon pairs unnatural for initiation were discriminated 
against by these methylations. For example, deficiency 
of RsmB and/or RsmD resulted in increased initiation 
by tRNA^uQ from CAC or CAU codons. Also, in our 
assays, where the native tRNAj^J was used to initiate 
from AUA codon, it initiated better in the absence of 
G966 methylation. We should mention that although the 
AUA:CAU pair with a CxA mismatch at the third 
position of the initiation codon is used as initiation 
codon in mitochondria, its use is facilitated by modifica- 
tion (formyl group) of the C34 of the CAU anticodon (3). 
In bacteria, where the anticodon of initiator tRNA is not 
modified, G966 methylation facing the CxA mismatch, 
may in fact be destabilizing. Interestingly, AUA 
decoding is special even at the step of elongation where 
C34 in tRNA^u is modified with agmatidine in archea 
(34) and lysidine in eubacteria (35). 

S9 tail contributes to the fidelity of initiation 

The S9 C-terminal tail suppressed initiation with 
tRNAg^ using CAC and CAU codons, the pairs 



Table 2. Summary of initiation activities (fold changes) of various 
anticodon stem mutants in the S9A3 strains expressed with respect to 



the parent E. coli SA 


Anticodon 


Initiation 


Anticodon stem 


Fold changes in 


in 


codon in 


mutation 


initiation in S9A3 


tRNA™ 61 


CAT 


in tRNA™ 61 


strain w.r.t. the 








control (SA strain) 






1GC (AU) mutant 


1.37 ± 0.3 






1GC (GU) mutant 


1.43 ± 0.3 


CUA 


UAG 


3GC mutant 


2.3 ± 0.13 






2GC (AU/GU) mutant 


3 ± 0.4 






No mutation in stem 


1 ± 0.19 



Table 1. Summary of initiation activities (fold changes) of various codomanticodon combinations in the S9A3, ArsmD and S9A3 ArsmD strains 
expressed with respect to (w.r.t.) the parent E. coli SA 

Initiation codon in CAT Fold changes in initiation in the given strain w.r.t. the control (SA strain) 



S9A3 ArsmD S9A3 ArsmD 



CAC 2.5 ±0.4 1.8 ±0.1 2.5 ± 0.05 

GUG CAU 3.8 ± 0.6 2.9 ± 0.3 4.5 ± 0.7 



GUA 


UAC 


1 ± 0.03 


1.52 ± 0.04 


2 ± 0.26 


GAU 


AUC 


1.1 ± 0.17 


0.85 ± 0.08 


0.96 ± 0.1 


CUA 


UAG 


1 ± 0.2 


1 ± 0.11 


1 ± 0.19 




AUG 


1.14 ± 0.08 


1.2 ± 0.09 


1.17 ± 0.06 




GUG 


0.6 ± 0.06 


1.2 ± 0.3 


0.67 ± 0.04 




UUG 


0.47 ± 0.05 


1.26 ± 0.04 


0.63 ± 0.1 




CUG 


0.76 ± 0.05 


1.4 ± 0.02 


0.9 ± 0.2 


CAU 


AUA 


1 ± 0.044 


2 ± 0.02 


1.35 ± 0.08 




AUU 


0.83 ± 0.1 


1.17 ± 0.25 


1 ± 0.08 




AUC 


0.98 ± 0.07 


1.2 ± 0.15 


0.95 ± 0.06 




ACG 


1 ± 0.13 


1 ± 0.02 


1.2 ± 0.12 
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unnatural for initiation. Similarly, initiation by tRNA^ e ^ 
was also prevented by the S9 tail. Essentially, in the 
absence of the S9 tail, G34 in place of C34 was 
accommodated better, which in turn led to better initi- 
ation. The S9 tail did not impact initiation with various 
other non-canonical codons (with mismatched first, 
second or the third positions of the codon) using the 
native initiator tRNA^jj. However, initiation with the 
naturally used alternate initiation codons (UUG, GUG 
and CUG) decreased, suggesting that S9 facilitated stabil- 
ity of codomanticodon interaction in these cases. Such a 
suggestion is supported by the MD simulations, which 
revealed that binding of tRNA™ 6 * could be subtly 
affected on deletion of the S9 C-terminal tail (Figure 7). 
Although such an effect is unlikely to make a significant 
contribution (as a fraction towards the overall affinity) to 
a fully matched wild-type codomanticodon pair 
(AUG:CAU), its contribution would mean more to a mis- 
matched pair (such as in the cases of naturally used alter- 
nate initiation codons). Second, the binding of tRNA^^ 
was shown to decrease only modestly on S9 tail deletion 
(8). The modest decreases in the binding of tRNA^ 61 may 
become significant for initiation from UUG, GUG and 
CUG. In this context, an unexpected observation has 
been that the absence of the S9 tail did not impact initi- 
ation (using tRNA^j) even from the codons mismatched 
at their third positions. 

Further, as revealed by the MD simulations (Figure 7, 
panels ii and iii), the regain of conformational dynamics 
and/or constrained hydrogen bond (between G:U pair) 
adequately explain increases in initiation from CAC and 
CAU codons by tRNAg^. The MD simulations also 
offered an explanation for negligible change in initiation 
with AUC:GAU pair (Figure 7, panel iv). Here, the 
H-bonds made between G34 of tRNA and the S9 tail 
are the same as reported for wild-type tRNA in the 
crystal structure (4). Thus, inspite of having G34 in its 
anticodon, this tRNA is largely insensitive to the 
presence or absence of the S9 tail. Finally, insights into 
the role of the S9 tail in decoding during initiation were 
gained by analysis of initiation with the 3GC variants 
of tRNA™ 61 where the mutants lacking either 2GC or 
all 3GC base pairs initiated significantly better in the 
absence of the S9 tail. This observation contrasted the 
earlier report (8), where the elongator tRNAs (tRNA phe , 
tRNA Cys ) having 2GC base pairs in the anticodon stem 
were preferentially stabilized in the absence of S9 tail. This 
may be explained by the fact that the in vitro study (8) did 
not use IF 3, whereas in the in vivo context, IF 3 would have 
a role to play in discriminating against the initiator tRNAs 
lacking the three GC pairs (36). 

Interestingly, the S9 tail phenocopies the role of IF3 in 
discriminating against tRNAs lacking the 3GC base pairs 
(36), suggesting that the S9 tail and IF3 may function in 
concert. The two key roles of IF3 are to discriminate start 
codons AUG, GUG and UUG from the rest and to dis- 
criminate initiator tRNA via its 3GC base pairs. On 
deletion of the S9 tail both the functions are affected, sug- 
gesting that the deletion of the S9 tail has effects in 
common with the deficiency of IF3; it would be of 
future interest to investigate these mechanisms further. 



P-site ribosomal proteins and the initiator tRNA 

This study highlights the functional/evolutionarily import- 
ant roles that the proteins located within the ribosomal 
P-site may have. It also rationalizes the selection and con- 
servation of the anticodon of the initiator tRNA. Initiator 
tRNA has a conserved C34 in its anticodon, which seems 
to be important for the ribosome to distinguish it from 
other tRNAs, allowing its selection from a pool of tRNAs 
richly populated by elongators. This process is facilitated 
by the P-site features, such as the methylated G966 and 
the S9 tail of the ribosome. The 3GC base pairs too are 
important in the process of selective targeting of the initi- 
ator tRNA to the P-site. However, the interactions 
between P-site features and the initiator tRNA seem 
complex and specific for various codon: anticodon 
interactions. 
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